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Improved Wilson QCD simulations at light quark masses 

A. C. Irving a , UKQCD Collaboration* 

a Theoretical Physics Division, Department of Mathematical Sciences, University of Liverpool, PO Box 
147, Liverpool L69 3BX, UK. 

We present preliminary results from UKQCD simulations at light quark masses using two flavours of non- 
pertubatively improved Wilson fermions. We report on the performance of the standard HMC algorithm at these 
quark masses where m^/m p < 0.5 in comparison with simulations using improved staggered quarks. 



1. SIMULATION PERFORMANCE 

Since 1997, UKQCD has been conducting dy- 
namical fermion simulations on a Cray T3E ma- 
chine using typically 128 processors (35 Gflops 
sustained). This has permitted production stud- 
ies on lattices up to 16 3 x 32 and m w /m p down to 
around 0.60 [Q. Study groups within the collab- 
oration are now planning physics programmes for 
the advent of the QCDOC multi-TerafLops ma- 
chine which is expected to be our main resource 
from 2003 onwards. Actions and algorithmic per- 
formance are key concerns. 

1.1. HMC performance 

The next figure shows estimates of the work 
(in Teraflops) required to obtain a 'new' configu- 
ration using three of the main methods currently 
available. 
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The data for the plot was mostly taken from 



that provided last year at LAT01 for the ses- 
sion on algorithm performance [g] . The top curve 
corresponds to the UKQCD 'clover- Wilson' esti- 
mate: 
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The accompanying curves correspond to Wilson 
fermions as measured by the SESAM collabora- 
tion j|] and the Asqtad MILC code ('order a 2 
tadpole' improved staggered fermions using the 
R- algorithm Q). The dashed curve corresponds 
to the Asqtad data but where the lattice spacing 
has been increased by 50%. The Asqtad scal- 
ing curve was taken from the performance data 
in Q. We performed an independent check by 
carrying out Asqtad simulations with 2 flavours 
on the same size lattice and machine used to pro- 
duce the UKQCD data. For this comparison, a 
'new configuration' amounts to 20 MD trajecto- 
ries for the Wilson and clover HMC runs and 6 
for the Asqtad runs. These numbers are based 
on using at least twice the measured integrated 
plaquette autocorrelation time in each case. The 
motivation for the extra comparison of Wilson at 
lattice spacing a with Asqtad at spacing 1.5a is 
the observation that the improvement for the lat- 
ter is much greater (0(a 2 a s )) than for Wilson at 
a given spacing a. Since the work scales like a~ 8 
this contributes very significantly to the saving in 
work at fixed physical lattice size and m^. 

The Wilson/SESAM simulation || appears 
to be 8 times faster than the improved Wil- 
son/UKQCD one. Detailed comparison shows 
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that a factor of 3.4 comes from solver iterations 
(better preconditioning and solver algorithm and 
better underlying condition number); a factor of 
1.8 from larger MD step size (higher acceptance, 
physics); a factor of 1.3 from the numerical over- 
head of the clover term. These naive comparisons 
do not take into account other factors such as rel- 
ative scaling behaviour and accessibility and effi- 
ciency of measurement operators. 

1.2. HMC stability 

These HMC runs are expected to be susceptible 
to instabilities (occasional large AH values and 
zero acceptance) when the fermion force term gets 
too large |6| . We have observed this effect directly 
on the clover- Wilson run at k = 0.1358 where we 
find it necessary to use a step size of 1 /400. When 
the step size is this, or smaller, we find that 64 bit 
for our field storage and matrix-vector manipula- 
tions is required so as to avoid a dramatic loss of 
acceptance due to rounding errors. This is so even 
though we always use full 64-bit arithmetic and 
sophisticated summing techniques in the global 
summing of our AH (energy difference) calcula- 
tions. 

1.3. Run parameters 

The main parameters of the new light quark 
run are shown in the following table. In physi- 
cal units, L « 1.5/m and m n « 430 Mev before 
extrapolation. 



trajectories/hour 
no. trajectories 
no. configs 
t ac (plaquette) 


(5. 2,. 1358) 
3 

2440 

122 

6.9(6) 


r /a 

a (using r = 0.49 fm) 

amps 

mps/mv 

mp S L 


5.3(5) 

0.0925(9) fm 
0.207(5) 
0.43(2) 
3.3(1) 



The integrated autocorrelation time for the pla- 
quette follows the previously observed trend that 
it decreases with decreasing quark mass. For 
comparison, r^=16(3) at k — .1350 where 
n r i'Ps/ m v ~ 0.7. Although not expected, this 
trend can be accommodated in simple models |l]] . 



The ratio l m n /m p ' is below the threshold for 
the continuum decay p — > irir but not for decay 
on this size of lattice. We are using maximum 
entropy methods to look for evidence of decay 
of the lattice state at the relevant energy. 

2. MEASUREMENTS 

2.1. Hadron masses and finite size effects 

Since m^L = 3.3(1) (significantly below 5) we 
anticipate finite size effects. For small and in- 
termediate size boxes there is evidence Q that 
hadron mass corrections follow a power law (£ -3 ) 
rather than exponential (exp (— L/Lq)) size de- 
pendence. A simple phenomenological model for 
this might be 



Sm h 



Ch 



{m^Lf 



(2) 



The dimensionless parameter Ch is presumably re- 
lated to the form factor of the hadron h. 

The figure below shows a compilation (from 
UKQCD and QCDSF §) of pseudoscalar and 
vector masses (fixed (3 = 5.2 and n s — n v ). The 
filled points are the raw masses extracted from 
multi-channel fits. 
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Aside from the point at very heavy quark mass 
(k = 0.1342), the mps data (filled circles) are 
well fitted by lowest order chiral PT (linear in 
1/k). 

We have estimated c,~6 and c p ps 7 from the 
finite size (FS) effects reported by the JLQCD 
Collaboration, who used the same action (clover- 
Wilson) |h]] on a variety of lattice sizes. We 
have checked that the corrections implied by these 
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values are consistent with the FS effects which 
UKQCD has previously measured using a slightly 
different action |T]j] . Using this model to 'cor- 
rect' the measured values, we get the open points 
shown in the figure. 

The FS-corrected pseudoscalar points are still 
well fitted by a straight line and yield K cr it — 
0.135978(9)(55). The first error is statistical and 
the second systematic - estimated by including 
the heaviest quark point in a quadratic fit. Taking 
the above FS corrections for the pseudoscalar and 
vector into account we estimate mpg/my w 0.45 
for the lightest quark simulation. 

2.2. Static potential 

The next figure shows the difference of the mea- 
sured static potential from the universal bosonic 
string model for all UKQCD data sets at fixed 
(3 = 5.2 |l| including this latest run (stars) at 
re = 0.1358 and a QCDSF set at re = 0.1340 {§. 
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Also shown is the corresponding result from a 
modest run we have made using the MILC Col- 
laboration's code at a « 0.130 fm [Q. Viola- 
tions of rotational symmetry are marked for r < 
0.7r w 0.35 fm for all these coarse lattices. Al- 
though the improved staggered fermion measure- 
ments (Asqtad) are at a coarser lattice spacing, 
they show smaller violations of rotational sym- 
metry than the clover- Wilson data. This is pre- 
sumably due to the significantly improved gluonic 
action. The Asqtad data (at mps/mv ~ 0.50) 
show a greater effect of the running coupling in 
the Coulomb term at short distances than the 
corresponding clover- Wilson data (lower at short 
distances). Within the fixed (3 clover- Wilson set 



there is no very clear trend with decreasing quark 
mass. The corresponding plot for clover- Wilson 
data at fixed lattice spacing |jj shows more clearly 
the effect of the running coupling. 

3. CONCLUSIONS 

We have analysed data from more than 2200 
trajectories of a clover-improved Wilson simula- 
tion at m q /m s 0.3. We have estimated finite 
volume corrections to meson masses. Light quark 
masses are accessible using clover-improved Wil- 
son fermions and standard HMC but simulations 
are slow and can suffer from instabilities. Im- 
proved staggered fermions appear to offer a more 
efficient route to configurations with light dynam- 
ical quarks, but physics measurements are in gen- 
eral more complicated. 
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